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A model consisting of random interacting anti-ferromagnetic (AF) grains coupled to a ferromag- 
netic (FM) layer is developed to study the exchange bias phenomenon. This simple model is able 
to describe several exchange bias behavior observed in real materials. Shifts in hysteresis loops are 
observed as a function of cooling field and average grain size. We establish a direct relationship be- 
tween cooling field dependence of exchange bias, coercivity and magnetization state on the AF-FM 
interface. We also verify that the exchange bias field is inversely proportional to the grain size, and 
this behavior is independent of the inter-grain interactions, AF/FM coupling and cooling field. 



When a ferromagnetic (FM) material is coupled to an 
anti-ferromagnetic (AF) material, under suitable condi- 
tions, a unidirectional anisotropy is observed. This re- 
sults in a shift in the hysteresis loop called exchange 
bias Because of its application to spin valves, ex- 

change bias has been studied extensively, but the roles of 
many parameters, such as magnetic domains and cool- 
ing field, have not been fully understood. There have 
been attempts to understand exchange bias in the mi- 
croscopic level [3- E|- In particular, recent publications 
on the domain state model [j| argued that exchange bias 
is due to different domain orientations in the AF bulk. 
In this paper, we make a theoretical analysis based on 
a microscopic model to have a better understanding on 
the roles of domain orientations. In contrast with pre- 
vious works 0, El, where domains were introduced by 
dilution, we model domains explicitly with grains on the 
AF materials. We derive a set of mean-field equations 
that consider a free boundary on the surface of each AF 
grain and an effective field on the surface due to interac- 
tions with other grains. We study the cooling field and 
grain size dependences of exchange bias. 

Our model for studying exchange bias consists of one 
FM layer on a square lattice coupled to multiple AF lay- 
ers on a cubic lattice. We used eight AF layers in all 
our calculations. Fig. ^ shows a schematic of the AF 
grains and defines the coordinate system which we shall 
use hereafter. Periodic boundary conditions are used for 
in-plane directions, and free boundary conditions are con- 
sidered in the z direction. The Hamiltonian for the FM 
layer is given by 

H PM = -J PM [d(s* f + hsf + Jzsfo-i] (1) 

(ij) 1 

where s*j are Heisenberg spins and J FM is the exchange 
coupling between FM spins. The first sum is performed 
over nearest neighbors (ij) within the FM layer. The i 
sum is taken over all lattice sites. The uniform external 
field h is chosen to be parallel to the easy axis x, and 
d = 0.1J FM is the anisotropy constant. J 1 is the exchange 
coupling between FM and AF layers, and several different 
values of J, are used in our calculations. For spins of AF 
layers we take Ising spins <Tj = ±1. The Hamiltonian for 
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FIG. 1: A schematic of FM layer with AF grains 
one AF grain is 

H main = - Jaf ^ ai(7 3 -^Y^Oi- K' a n (2) 
(ij) i n 

where the coupling J AF is chosen as — 0.5J FM - The first 
and second sums are taken over the nearest neighbors 
and all lattice sites of each grain, respectively. /i^ ur is an 
additional effective field on the surface, and the n sum 
is performed over the surface of the grain. The effective 
field h^ r consists of two terms. The first term is the con- 
tribution from the nearest neighbor spins belonging to 
different grains, J s J2j 17 ji where we take J g = — 0.15J FM 
(| JJ < | J A p I)- The second term is J z sf if the site n be- 
longs to the interface layer, and k is the nearest neighbor 
site of n in the FM layer. 

To generate the grain configuration, we have take the 
following procedure. We start from an "unoccupied" AF 
lattice and insert n s seeds, where each seed is the starting 
point of a distinct grain. From these seeds, we grow 
grains by the following process. Randomly choose a site. 
If an empty site is chosen, choose another site. If the 
chosen site belongs to a grain, let all empty neighbors 
of this site belong to the same grain. Repeat this step 
until all sites are occupied. Fig. [21 shows the grain size 
distribution with 512 grains on 32 x 64 x 8 lattice averaged 
over 512 samples. Grain configurations with 512 grains 
and average grain size of 32 lattice sites are used unless 
otherwise stated. 
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FIG. 2: Plot of grain size distribution with mean grain size = FIG. 3: Hysteresis loop for Ji = —0.5 (white squares) and 
32. J 1 = 0.5 (black circles), with representative error bars ob- 

tained from 64 independent simulations. 



We derive a set of coupled mean-field equations for 
the magnetization of each AF grain. These equations 
account for the free surfaces and the surface fields due to 
interactions with other grains. Details of derivations will 
be shown elsewhere 0. For the jth grain, it follows that 



For the FM layer the free energy is calculated as 
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where M is the FM magnetization, q FM — 4 is the coor- 
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dination number of FM layer, h k 

defined by hf* = {h + J : m^ /x )e x 
integral to be evaluated numerically, 



is an effective field 
. I[M,ht /x ] is an 



Psjkls tanh[/3( J AF q s m^ + h + J s q s mt)} (3) I[M, h] = ds exp (3(J FM q FM M + h)-s + f3d{s x ) 
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where the variables ip and x are used to represent sublat- 
tices, that is, rrij and are the magnetization of one 
sublattice of the jth grain and that of the other sublat- 
tice, respectively. For each grain, there are a pair of sym- 
metric equations for both sublattices. As a total, Eq. © 
consists of a set of 2n g equations. The variables in Eq. © 
are defined as follows: q = 6 is the coordination number 
in the bulk, q s is the average coordination number on the 
surface, q g is the average coordination number between 
sites belonging to different grains, and (3 is the inverse 
temperature. By defining the densities of bulk/surface 
sites as the number of bulk/surface sites divided by the 
total number of sites in the simulation box, we express 

the mean-field equation in terms of densities. and 
ptj X are tne densities of bulk sites and surface sites for 
sublattice ip/x> respectively. Pg !lt , g j is the density of free 
surface, P5 FMJ - is the density of sites adjacent to the FM 

layer, and Pgj k is the density of surface sites adjacent to 
the neighboring fcth grain. The first term in the right 
hand side of Eq. (J3J) gives the contribution to m* due to 
magnetization ordering in the bulk. The second term 
gives the contribution due to the free surfaces at z = 0. 
The third term is due to interactions with the top FM 
layer and the last term is due to inter-grain interactions. 



(5) 

The x component of M, M x , appears in Eq. J3J). Com- 
putation of mean-field magnetizations consists of solv- 
ing Eq. @ iteratively and simultaneously minimizing 
Eq. Q. 

Let us explain the field-cool process. From now on, we 
use dimcnsionless units with J FM = 1, and the tempera- 
ture T is given in units of J FM /k B , where k B is the Boltz- 
mann constant. The system is first cooled under the field 
h coal from T = 5 to T = 0.2 in steps of T l+1 = 0.95Tj. 
There are two possible domain orientations for the AF 
grains, with sublattice ip pointing in the +1 direction 
and x pointing in the —1 direction, or vice versa. Dur- 
ing the field-cool process, we attempt to obtain the global 
free energy minimum by considering two different domain 
orientations for each AF grain and pick up the one with 
lower free energy. At the end of the field-cool process, the 
temperature is kept constant and the field h is reduced in 
steps of 0.01 and with a higher resolution near the region 
of magnetization reversal (0.002). At each step, the mag- 
netization is calculated using the previous magnetization 
state as the starting values for Eq. and Eq. In this 
way, meta-stable states can be traced out to obtain the 
hysteresis loop. We simulate the field-cool process and 
the hysteresis loop based on the mean-field equations. 

Fig-Glshows the hysteresis loops at T = 0.2. Hysteresis 
loops shifted to the right for J t = —0.5 (white squares) 
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FIG. 4: The dependence of H EB and H c on cooling field for 
Ji = —0.5 (white squares) and Ji = 0.5 (black circles). The 
bottom plot shows the magnetization m and staggered mag- 
netization m* at the AF/FM interface at T = 0.2. 
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FIG. 5: Grain size dependence of exchange bias for Jj = 
—0.5, J g — —0.15, ftcooi = 1.0 (filled squares), Ji = —0.5, J g = 
— 0.15, /icooi = 2.0 (empty circles), Ji = — 0.5, J g — 0,h caol 
I . I ! (filled diamonds). Error bars are smaller than the size of 
the symbols. 



and left for J 1 = 0.5 (black circles). Field-cooling is per- 
formed at h cool = 2. The origin of the shift comes from 
excess magnetization in the AF bulk created during field- 
cooling. This result is consistent with the intuitive pic- 
ture discussed by Nogues and Schuller 0. The shift in 
hysteresis loop should depend on the magnetization state 
immediately before the hysteresis loop is traced out. We 
shall modify the magnetization state by using different 
cooling fields. 

We performed a systematic study on the dependence 
of the exchange bias field H EB and coercivity H c in 
the field-cool process. H EB is defined as (H+ + H_)/2 



where H + and i/_ are the fields at which the magneti- 
zation is zero in the hysteresis loop. Similarly, we define 
H c = (iJ + — H-)/2. Fig. ^] shows that H EB increases 
from zero to a maximum at about /i cool = 2.5 and de- 
creases at hcool > 6. The hysteresis loop remains un- 
shifted with ft, cool = when cooled from a demagnetized 
state; this observation is consistent with experiments on 
exchange bias pj, |S|, |9| . Increase and decrease of H EB and 
H c with cooling field is also reported in real materials 
such as Ni/NiFe 2 4 [13 and permalloy/CoO [Hill. As 
shown in the bottom plot of Fig. 01 at h < 2.5, the spins 
of AF-FM interface are locked in a magnetization state 
with staggered magnetization m* rs 1, and at h > 6, 
the staggered magnetization vanishes. Our result shows 
that cooling field of /i cool = 2.5, for the sets of parame- 
ters used in the present study, produces maximum excess 
magnetization in the AF-FM interface while maintaining 
non-zero staggered magnetization, causing a maximum 
exchange bias. At higher cooling field, spins in the AF 
material are aligned during field-cooling; when the field 
is decreased to sweep the hysteresis loop, the AF grains 
get locked in AF states with random domain orientations 
that generate little excess magnetization. Hence there is 
a close relationship between the magnetization state of 
the AF/FM interface and H EB . We also found, as in [Tlj . 
that there is a decrease of coercivity corresponding to an 
increase of exchange bias. The dependence of exchange 
bias on cooling field is of particular importance in fabri- 
cation of exchange bias devices where tuning of exchange 
bias is desirable. 

Fig. [SI shows grain size dependence of exchange bias. 
Average grain sizes of 8, 16, 32, 64 and 128 are used. Sim- 
ulations with interacting/non-interacting grains and dif- 
ferent cooling fields were performed. The exchange bias 
field is inversely proportional to the grain diameter D, 
H EB ~ l/D. For non-interacting grains (filled diamonds) 
and strong cooling field (empty circles) , plots of exchange 
bias versus inverse grain diameter D^ 1 fall on a straight 
line very well. For smaller cooling fields (filled squares), 
exchange bias goes to zero asymptotically at large grain 
sizes. Simulations with Ji = —0.75 were also performed 
and the exchange bias field changes only slightly com- 
pared to simulations with J/ = —0.5. This inverse rela- 
tionship of exchange bias on grains size is also reported 
in exchange bias systems of permalloy/CoO bilayers [T^ 
and Cr 7 oAl 3 o/Fei9Ni 8 i bilayers Q. 

Other models has been developed to describe the 
mechanism of unidirectional anisotropy. Stiles and 
McMichael O used an ordered granular model to ex- 
plain exchange bias through partial domain wall forma- 
tion. Scholten et al. |6J used mean-field equations on local 
magnetization to explain the domain state model. The 
model presented here is distinct from all previous models. 
The mean-field equations for an explicit grain distribu- 
tion has been given for the first time in this paper, and 
the domain state model is realized explicitly by actual 
grain distribution. 

To summarize, we developed a simple model that cap- 
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tures many features of real exchange bias systems. We 
found a direct relationship between the cooling field de- 
pendence of the exchange bias, coercivity and magneti- 
zation states on the AF-FM interface. We also verified 
that the exchange bias field is inversely proportional to 
the AF grain sizes and this relationship is independent 
of the inter-grain interactions, AF/FM interactions and 
cooling fields. Lastly, we would like to mention that the 
simulations based on mean-field equations used in this 



paper are general, and may be used to study other ex- 
change bias parameters, such as surface roughness, per- 
pendicular coupling and film thickness. 
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